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ABSTRACT 

We study the scattering of low-energy Cosmic Rays (CRs) in a turbulent, compres- 
sive MHD fluid. We show that compressible MHD modes - fast or slow waves with 
wave lengths smaller than CR mean free paths induce cyclotron instability in CRs. 
The instability feeds the new small-scale Alfvenic wave component with wave vectors 
mostly along magnetic field, which is not a part of the MHD turbulence cascade. This 
new component gives feedback on the instability through decreasing the CR mean 
free path. We show that the ambient turbulence fully suppresses the instability at 
large scales, while wave steepening constrains the amplitude of the waves at small 
scales. We provide the energy spectrum of the plane-parallel Alfvenic component and 
calculate mean free paths of CRs as a function of their energy. We find that for the 
typical parameters of turbulence in the interstellar medium and in the intercluster 
medium the new Alfvenic component provides the scattering of the low energy CRs 
that exceeds the direct resonance scattering by MHD modes. This solves the problem 
of insufficient scattering of low-energy CRs in the turbulent interstellar or intracluster 
medium that was reported in the literature. 
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1 INTRODUCTION 

] Cosmic rays (CRs) and magnetic fields are essential compo- 

■ nents for many astrophysical ecosystems, including galaxies 

■ and clusters of galaxies (see Schlickeiser 2002). In many in- 
stances, e.g. Milky Way, the pressure of CRs and magnetic 
fields is larger than the gas pressure. As a rule, astrophysical 
magnetic fields are frozen in turbulent plasma and move to- 
gether with it. As a result, CRs interacting with turbulent 
magnetic fields get scattered and accelerated (see Melrose 
1968, Schlickeiser 2002). 

The magnetohydrodynamic (MHD) approximation is 
widely used to describe the actual magnetized plasma tur- 
bulence over scales that are much larger than both the mean 
free path of the particles and their Larmor radius (see Kul- 
srud 2004). The theory of MHD turbulence has become 
testable recently due to numerical simulations (see Biskamp 
2003) and this provided reliable foundations for describing 
turbulence-CRs interactions. The simulations (see Cho & 
Lazarian 2005 and ref. therein) confirmed the prediction of 
magnetized Alfvenic eddies being elongated along magnetic 
field (see Shebalin, Matthaeus & Montgomery 1983, Higdon 
1984) and provided results consistent with the quantitative 
relations for the degree of eddy elongation obtained in Gol- 
dreich & Sridhar (1995, henceforth GS95). 

Scattering of CRs is an essential part of both CR propa- 
gation modes and models of CR acceleration. Eflicient scat- 
tering is usually postulated (see Schlickeiser 2002), which 



ensures high degrees of coupling of CRs and magnetized 
plasma. In addition, eflicient scattering provides apprecia- 
ble second order Fermi acceleration and enables the return 
of CRs into the shock to ensure the first order Fermi accel- 
eration. 

This corner stone of CR physics has been challenged 
recently when it became clear that Alfvenic eddies are 
stretched along magnetic field direction. As the interaction 
between CRs and such elongated eddies is weak (see discus- 
sion in Lerche & Schlickeiser 2001), this resulted in the pre- 
diction of long mean free paths for Milky Way CRs (Chan- 
dran 2000, Yan & Lazarian 2002, henceforth YL02). YL02 
and Yan & Lazarian (2004) attempted to remedy the sit- 
uation by appealing to CR scattering by isotropic sound- 
like fast modes. However, plasma-dependent damping of fast 
modes made the scattering very different in different parts 
of the interstellar medium. Is such a radical change of the 
CR scattering picture absolutely necessary? 

We note, that the problem of CR scattering goes well 
beyond the Milky Way physics. Brunetti (2006) discussed 
the implications of suppressed CR scattering on the acceler- 
ation of CRs in the clusters of galaxies. Is there any process 
through which scattering by fast modes can provide high 
efficiency of CR scattering? Below we consider such a pro- 
cess that is related to CR feedback on MHD turbulence. 
We show that compression of CRs induces instability that 
results in the generation of modes that are parallel to the 
magnetic field. Such modes that are also frequently referred 
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to as slab modes have been long employed in the models of 
CR propagation (sec, e.g., Jokipii 1966), their origin, how- 
ever, was somewhat mysterious. This paper provides a phys- 
ically motivated mechanism for the generation of slab modes 
and quantify the efficiency of their generation. 

In what follows we discuss the properties of compress- 
ible MHD turbulence in § 2. We describe the kinetic in- 
stability that develops in CRs when the magnetic field is 
compressed on scales less than the CR mean free path in 
§ 3. We consider the non-linear saturation of the instability 
in § 4 and its large-scale cut-off that follows from the inter- 
action of the instability waves with the ambient turbulence 
in § 5. The feedback of CRs to compressible turbulence is 
considered in § 6. The implications of our work for various 
ISM phases and intrarcluster medium (ICM) are considered 
in § 7. The short summary of this work is presented in § 8. 



2 COMPRESSIBLE MHD TURBULENCE 

In this section we provide a summary of the current knowl- 
edge of compressible MHD turbulence that we appeal to in 
our work. As we mentioned earlier, we describe an alterna- 
tive mechanism for the scattering of low-energy CRs that 
provides rather small mean free paths for CRs. At scales 
larger than the mean free path CRs arc well coupled to the 
flow of the thermal plasma. Therefore to identify new effects 
we concentrate on the MHD fluctuations on scales equal or 
smaller than the mean free path^. 

The GS95 model quantifies turbulence anisotropy, in- 
troducing the so-called critical balance relation fcy ~ k^/^ , 
where A;|| and k± correspond to, respectively, the parallel and 
perpendicular wavenumber of eddies measured in respect to 
the local magnetic field. This scaling is true for both Alfvcnic 
and pseudo-Alfvenic motions, which are the incompressible 
limit of the slow modes. GS95 and later Lithwick & Goldre- 
ich (2001) argued that the slow modes are passively advected 
by the Alfven modes, while the energy exchange between the 
modes is small. Numerical studios (Cho & Vishniac 2000, 
Maron & Goldreich 2001, Cho, Lazarian & Vishniac 2002, 
Miiller, Biskamp & Grappin 2003) provided reasonable sup- 
port for the critical balance condition^. We introduce the 
outer scale La and use GS95 scaling in the form 

fell ~ kf-'L-^^'\ 6v - VA{k^LA)'^\ (1) 

We feel that GS95 provides a good starting point for 
studies of mildly compressible, and oven supersonic MHD 
turbulence. Indeed, numerical calculations in Cho & Lazar- 
ian (2002, 2003, henceforth CL02, CL03) showed that scal- 
ings of the slow and Alfven modes in compressible MHD 
are similar to their scalings in the incompressible case. The 
fast mode perturbations, on the other hand, are found to 

^ A notable exception from this rule, MHD shocks, relevant for 
CR acceleration, is not considered in this paper. 
^ Note that the GS95 prediction of the Kolmogorov spectrum 
—5/3 spectral index stayed more controversial, as Maron Sz Gol- 
dreich (2001) reported the spectrum closer to —3/2. While the 
ongoing work (sec Muller & Grappin 2005, Boldyrev 2005, 2006, 
Beresnyak & Lazarian 2006) attempts to improve our understand- 
ing of Alfvenic turbulence, we shall use the original GS95 scalings 
for the sake of simplicity of our presentation. 



be mostly isotropic with a power-law index of about —3/2 
(see CL02), which is the index of so-called acoustic wave 
turbulence. The coupling of the fast and Alfven modes was 
shown to be weak, which allows separate studies of the cor- 
responding cascades provided that the Alfvenic turbulence 
is strong, i.e. it evolved to develop the critical balance. 

In the following treatment we will be primarily inter- 
ested in two manifestations of compressible MHD turbu- 
lence, one of which is perturbations of the magnitude of the 
magnetic field and the other is wave damping through cas- 
cading by the ambient Alfvenic turbulence. Let us briefly 
discuss how much of our results depend on the adopted 
model of MHD turbulence. 

As we are dealing with very small perturbations down 
in the inertial range, Alfvenic mode, having magnetic field 
perturbations perpendicular to the local magnetic field, has 
very little effect on the magnetic field magnitude, therefore, 
as far as the magnetic field compression is concerned, we are 
dealing with compressible modes. We introduce the power- 
law energy spectrum for the velocity perturbations of the 
compressible modes, E{k) ~ fc^", whore P = 5/3 for the 
Kolmogorov- type, and 3/2 for the acoustic type spectrum. 
Such a spectrum translates into velocities at the scale of I, 
which is 5vi ~ va{1/L)^ , where //, = {(3 — l)/2 and L is not 
necessarily the injection scale for MHD turbulence, but a 
quantity which contains information on both the injection 
scale and the efficacy of driving. This scaling might not be 
valid up to scales as large as L. If both sonic and Alfvenic 
Mach numbers are around unity at the injection, we expect 
L to be of the order of the actual injection scale. Efficacy 
of driving depend on its type, such as mostly solenoidal, 
or mostly compressive (supernova shocks, etc.). A precise 
estimate of the parameter L requires either detailed knowl- 
edge of the physics of driving, or direct measurement of the 
compressive fluctuation intensity somewhere in the inertial 
range. 

The Alfven wave damping by Alfvenic turbulence we 
use in § 4 assumes GS95 scaling and anisotropy with some 

outer scale La- This scale might correspond to the isotropic 
injection of energy at scale L and the injection velocity of va, 
i.e. the Alfven Mach number Ma = {5v/va) = 1- This model 
can be easily generalized for both Ma > 1 and Ma < 1 at 
the injection. Indeed, if Ma > 1, instead of the driving scale 
Ld for La one can use the scale at which the turbulent ve- 
locity gets equal to va- For Ma ^ 1 magnetic fields are not 
dynamically important at large scales and the turbulence 
follows the Kolmogorov cascade vi ~ l^^^ over the range of 
scales [Ld,La]- This provides La ~ LdM^'K If Ma < 1, the 
turbulence obeys GS95 scaling (also called "strong" MHD 
turbulence) not from the scale Ld, but from a smaller scale 
I' ~ LdMa (Lazarian & Vishniac 1999), while in the range 
[Ld, I'] the turbulence is "weak". The velocity at scale I' is 
expressed as U;/ ~ vaM'a, so that the "effective" value of 
La wiU be La = LoMj". 

All in all, given the strength and the nature of driving 
in a particular astrophysical environment one may estimate 
the two parameters, L and La that determine the veloc- 
ity perturbations of the compressible and Alfvenic mode at 
small scales. 

In § 3 we deal with the fluctuations of the magnetic field 
magnitude squared. The normalized amplitude of these fluc- 
tuations denoted as A will depend on the plasma /J which is 
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the ratio of the gas pressure to the magnetic field pressure. In 
high-/? plasmas will be perturbed mostly by slow waves 
(CL03) and the value of A will be equal to 2{5v/va) sinO, 
where 6 is the angle between the wave vector and the mag- 
netic field. In the inertial range of strong turbulence the slow 
mode exhibits the same anisotropy as the Alfvcnie mode. 
Therefore 6 is close to 90°, so we can disregard the angular 
factor, i.e. A = 25v/va- The perturbation made by the fast 
mode in high beta plasmas is smaller by a factor of va/cs, 
where Cs is sound velocity. In low-/? plasmas this situation 
is reversed, with the slow mode only marginally perturbing 
B, but we may use the same expression A — 2{5v/vA)sm9 
for the fast mode. We estimate the angular factor as of the 
order of unity, since fast modes are almost isotropic (CL02). 
In other words, the expression A = 25v/va is interpreted 
as the compression factor for the most compressive mode 
which is the slow wave in high-/? plasmas and the fast wave 
in low-/? plasmas. 

The lower limits to the scales we described are deter- 
mined by the damping of MHD modes. While the Alfvenic 
mode in fully ionized media is supposed to be damped 
at the thermal Larmor radius, the compressive modes are 
damped by more efficient coUisional and coUisionless damp- 
ing (Ginzburg 1961, Barnes & Scargle 1973, YL04). In this 
paper we introduce the value of the most compressive mode 
cut-off scale as Zcut- 



3 INSTABILITY OF COMPRESSED CRS 

It is obvious that at scales less than their mean free path CRs 
can be treated as a coUisionless fluid. Particles in the coUi- 
sionless fluid preserve the adiabatic invariant p\/ B, where 
B is the the magnetic field strength and p± is the momen- 
tum perpendicular to the magnetic field. MHD compressive 
modes change the magnitude of B so that the distribution 
in momentum space becomes anisotropic (see Chew, Gold- 
berger & Low, 1956). 

Such a distribution is subject to a number of instabili- 
ties, some of which are hydrodynamic, i.e. involve the change 
of the entire distribution function, while others are kinetic, 
i.e. involve a chance in a fraction of particles, that is reso- 
nant with a particular wave-mode. Well-known examples of 
hydrodynamic instabilities are firehose and mirror instabil- 
ities (see e.g. Mikhailovskii, 1975). Hydrodynamic instabili- 
ties are typically fast with the largest wavenumber growing 
almost as fast as the gyrofrequency, but have a threshold, 
i.e. small deviations from isotropy do not induce instability. 

While compressive motions can generally induce rather 
large changes in the magnitude of B on scales at the injec- 
tion scale of turbulence L, in §3 we estimate mean free path 
and show that it is much smaller than L. Therefore the com- 
pressions of the magnetic field we deal with are too small to 
induce hydrodynamic instabilities. 

It has been well known that the momentum distribution 
functions with p^ = Py > Pz, are subject to kinetic instabil- 
ity called gyroresonancc instability (Sagdeev & Shafranov, 
1961, Mikhailovskii, 1975, Gary, 1993). This instability re- 
ceived less attention than its hydrodynamic counterparts. 
However, it is pretty fast for a power-law distribution of 
CRs, as we demonstrate below. 

For a power-law distribution of CRs the growth rate of 



the cosmic-ray- Alfven gyroresonancc instability (henceforth 
GI) can be estimated as (see Appendix): 



. , ncR(p> mws/fc||) 
7cr(k||) = ±iJpi AQ, 



(2) 



where ncR(p > mf2/fc||) is the number density of CRs with 
momentum larger than the minimal resonant momentum for 
a wave vector value of fcy , rn is the proton mass, n is the den- 
sity of the thermal plasma, tOpi is the ion plasma frequency. 
Q is a numerical factor, defined in the Appendix. The ± 
sign corresponds to the two MHD modes. We shall concen- 
trate on the Alfven mode, corresponding to the plus sign, as 
those are less subjected to linear damping (see § 2, § 6). As 
we will demonstrate in the next chapter, when anisotropy 
is created by compressive turbulence, the anisotropy factor 
A — {p± — p||)/p|| will be small and will change its sign on 
the scale of the mean free path, depending on two competi- 
tive mechanisms scattering which tends to isotropize mo- 
mentum distribution, and magnetic field compression which 
tends to make it oblate or prolate. 

We assumed that the unperturbed distribution of CRs 
is isotropic and follows a power law i.e. Fq ~ p^"^^ where 
a is conveniently defined as the power-law index for a one- 
dimensional distribution (or particle density). For example, 
around the Earth a ~ 2.6 up to the energies of 10^* eV. 
Note, that in order for the total energy to converge at high 
energies, a should be larger than 2. 

The expression for the instability rate, assuming A = 
25v/va (see § 2), could be written as 



Sv f T \ 



(3) 



where is a Larmor radius of a CR resonant with a partic- 
ular wave vector fcy = mQ,/p, ro is the 1 GeV proton Larmor 
radius and 



1 ■ 10 cm 
ncn{rp > ro) 



PC. (4) 



4 NON-LINEAR SUPPRESSION AND 
SATURATION 

We introduce the CR mean free path A below which CRs 

could be treated as coUisionless and the instability described 
in § 3 is active. In the absence of other scattering processes 
the CRs are scattered by the slab-type motions generated 
by the instability above. Let us estimate A following Longair 
(1994). If the change of magnetic field direction is </) ~ SB/B 
the scattering that is a random walk requires N ~ ^/(p^ 
interaction and 



(5) 



where we designated SB as the magnetic field perturbation 
pertaining to a particular wavenumber, i.e. SB'^ w E{k)k. 
We can consider 5B as a function of either k or the reso- 
nant Larmor radius (see Longair 1994). As the instability 
grows, 5B will grow, which reduces A. On the other hand, it 
is the mean free path A which determines the scale at which 
compressions of the magnetic field are important. This can 
be understood as follows: the CR distribution "remembers" 
the perturbed value of the magnetic field and its anisotropy 
only during the time the typical particle travels its mean free 
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7CR 



lAl) \ b ^ 



(6) 



path. Once particles scatter significantly, the anisotropy of 
the distribution is efi'ectively "reset". As a result only low 
amplitude motions on scales less than A excite the instabil- 
ity, or in other words, the degree of anisotropy A is deter- 
mined by the local perturbation of the magnetic field on the 
the scale A. We call this process non-lmear suppression. 

The instability grows as d(SB^)/dt = 7cr('5B^) where 
the injection of energy is happening at the scale of the mean 
free path, i.e. 

where eqs. Q and were used. We see that according 
to the above equations SB perturbations will grow as i^/^^ 
thus reducing A virtually to rp. In other words, the non- 
linear suppression is not able to constrain the development 
of instability and we have to consider other non-linear pro- 
cesses, such as wave steepening. 

Steepening does not occur for a monochromatic circu- 
larly polarized Alfven wave as the amplitude of the magnetic 
field stays the same. However, for a collection of waves with 
different wavelengths the amplitude of the magnetic field 
fluctuates in space and time and therefore the steepening 
effect is present. The steepening rate can be estimated as 



'Ystccp 



-{SB/Bfk^^VA, 



(7) 



where the "-" sign reflects the fact, that steepening damps 
the instability. 

By comparing and Q we get the equilibrium or 
saturated amplitude of the instability-induced perturbations 



B 



1/2 



^l/(2M+2)^M/(2n-t-2) 



(|J.-a+2)/{2^^+2) 



(8) 



which for a = 2.6 and /i = 1/3 produces a rather shallow 
spectrum of perturbations, E{k) ~ (SB)^ /k ~ fc"" *. 

Combining Eqs. JSJ and (|HJ| one gets that the energy of 
the slab modes at fc|| ~ Vp is supplied from the compressions 
at scale 



A ; 



£^l/(f'+l)^M/(M+ 



(a-l)/(^ + l) 



(9) 



So far we assumed that the turbulent compressible mo- 
tions are not damped. This is a good approximation until A 
is larger than the compressive mode cutoff scale Zcut. If, on 
the other hand, Zcut > A the compression for the instability 
is supplied from the eddies at the damping scale, namely, 
6v/va ~ {lcut/Ly^^{X/lcut)- The modification of our formu- 
lae is self-evident. Instead of eq. ||5J one gets 

1/4 



55 



1/2 



^l/4^A./4^(l-n)/4 



Vro / 



(3-a)/4 



(10) 



which, for the same value of a = 2.6, corresponds to a 
steeper spectrum of E{k) ~ k~^'^ . 



5 DAMPING BY ALFVENIC TURBULENCE 

The instability we considered in § 3 has the largest growth 
rates for the wave vector parallel to the field. This is due to 
the fact, that the phase of a resonance can be kept constant 
for a long time only if k± is small. The ambient turbulence 
non-linearly damps the instability through a process that 
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Figure 1. Decorrelation of a plane, kj_ = Alfven wave by tur- 
bulence. Lower picture shows the energy density of a wave in 
cylindrical k-space. In this picture Alfven waves were injected at 
= 17. Wave energy is being transferred in the direction of 
k±_ axis, which is typical for decorrelation by MHD turbulence. 
Wave decay was exponential after driving was switched off. Upper 
plot shows decay rate of the wave vs its wavenumber. A detailed 
discussion is presented in Beresnyak & Lazarian (2006). 



is analogous to the suppression of the streaming instability 
(YL02, Farmer & Goldreich 2004). In what follows we find 
the lower limit on k± using the approach similar to that in 
Farmer & Goldreich (2004). We also provide the results of 
numerical calculations that validate this approach. 

For Alfvenic turbulence we adopt the GS95 scaling, Q, 
which refiects the tendency of eddies to get elongated along 
the magnetic field. For the sake of simplicity we take the 
scale La, introduced in §2, equal to the scale L. This is not 
necessarily true for any astrophysical environment, however 
our formulae are trivially generalized for the case of La L. 

Consider a wavepacket of Alfven waves that moves 
nearly parallel to the magnetic field with the dispersion of 
angles (5fej_/fc|| ~ Ok- The individual waves follow the local 
direction of the magnetic field lines. As a result, the dis- 
persion in angles of the wave packet cannot be less than 
the dispersion of angles due to the ambient Alfvenic tur- 
bulence. Ok > 0bk- The latter for the GS95 model (see 
eq. lEJ) is Obk ~ SB^/Bq ~ (k^L)''^^^. The modes with 
minimal Ok are the fastest growing ones. As we establish 
below (see Eq. llU 'l. they are the least damped. Therefore 
for our simplified treatment we shall limit our attention to 



the wavepackets with resonant fc,. 



r„ and Obk ~ Ok- One 



can determine the characteristic perpendicular wavenumber 
kj_ ^ Sk± ^ rp^{rp/ LY^^ of the "most parallel modes" that 
are created by streaming CRs. 
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The strong Alfvenic turbulence decorrelates the wave- 
packet with k± on the time scale of v±k_i_. Thus using the 
above expression for k± and Eq. Q we get 

, ,2/3r-l/3 -l/2r-l/2 

7turb ~ "k±V± ~ ~VAkl L ' ~ -VATj, ' L ' , (11) 

which, up to the "-" that we used to denote the damping 
nature of the process, coincides with the damping rate ob- 
tained in Farmer & Goldreich (2004) and with the results of 
our numerical simulations shown in Fig. 1. 

A comparison between eqs. (|3J and IIH indicates that 
for the spectral index of CRs a > 3/2 the ambient Alfvenic 
turbulence provides an upper limit on the scale of perturba- 
tions that arise from compressible-induced instabilities even 
without accounting for nonlinear suppression. 

If we use nonlinear suppression, the critical scale can be 
obtained by using eqs. Q, (|HJ and lllll . For a > 5/3 our 
instability is damped for all scales, larger than 

rp.crit ~ ro [L '^r[, ) . (12) 

Therefore the spectrum of plane Alfven waves given by 
eq. © will protrude from rp^crit down to rp^min which cor- 
responds to minimum energies of CRs. 



6 FEEDBACK ON COMPRESSIBLE 
TURBULENCE 

Linear instability theory does not describe the energy trans- 
fer between particles and CRs. Moreover, in calculating lin- 
ear response we have not included the CR input into the real 
part of the dielectric tensor. Speaking of our description of 
CR scattering by waves, we can consider each resonant scat- 
tering event in the rest frame of the wave where it is a pure 
pitch-angle scattering, i.e. there is no energy transferred be- 
tween the particle and the wave. Basically, the model of 
linear instability and its nonlinear suppression and satura- 
tion developed in §4 is unable to fully describe the transfer 
of energy from compressible modes to CRs and from CRs to 
high-frequency waves which are finally dissipated by steep- 
ening. 

We would like to note that the transfer of energy is, in 
principle, might be very important for the particular mod- 
els of turbulence we are using in §2. They are so-called 
Kolmogorov-type models ^ in which the turbulence is ini- 
tiated at large scales, is being local in k-space and is char- 
acterized by the flow of energy from large to small scales*. 
Both slow wave passive transfer by strong Alfvenic turbu- 
lence and weak fast wave turbulence models we use in § 2 
are of this type. When energy is being drained out of such 
turbulence it does not any longer follow universal power-law 
scalings we described in §2. In this section we will estimate 
the transfer of energy from compressible motions and how it 
modifies the spectra of turbulence and processes described 
in §4. 

Let us first figure out the processes of transfer of energy 

^ This very general concept should not be confused with a partic- 
ular Kolmogorov theory for incompressible strong Navier-Stokes 
turbulence 

* In some cases this is a flow from small to large scales, as in 
Lcngmuir wave turbulence (see, e.g., Kaplan, Tsytovich, 1973). 



from CRs to waves. As we noted before, there is no transfer 
of energy in the frame moving with the wave. Therefore, if 
we have waves moving in only one direction, particles will 
pitch-scatter and establish the drift equilibrium with such 
waves. The drift velocity will be equal to the Alfven velocity. 
This is the case of the well-known streaming instability if we 
can neglect the wave damping. In our case, however, we have 
waves moving in both directions, so the particle can lose or 
gain energy in the lab frame by multiple scattering. If each 
scattering event changes pitch angle by <j!> ~ SB/B (c.f. eq 
(5)) the energy is changed by 2p±VaSB / B . This occurs at 
a Larmor frequency, i.e. c/vp. Now, starting with two above 
equations we can talk about two processes, one of which is 
the well-known diffusion in energy space, described by the 
diffusion coefficient Dpp estimated as 

This process is mostly describe the particle-wave equi- 
librium state where energy is slowly redistributed between 
particles. Suppose, however, that we efficiently drain en- 
ergy from the waves so that the situation become non- 
equilibrium. In this case we talk about the second pro- 
cess which is a directed transfer of energy from particles 
to waves. The rate of this process depends on the degree of 
non-equilibrium, but generally cannot exceed the inverse of 
the minimum time at which particles can lose the energy 
available from anisotropy, i.e. 

1 / VA SB I 

7cx = l/rox,miii ~ — (14) 

This rate is actually B/{5BA) times higher than a 
steepening rate. More careful estimate show that the vol- 
umetric rate of the energy exchange between particles and 
waves divided by the volumetric steepening energy rate 
is equal to B/{2SBAQ) which is much larger then unity. 
Therefore we can conclude that if the energy is drained from 
waves at a steepening rate, it will be efficiently resupplied 
by CRs. 

The anisotropy of CR distribution is, in turn, supplied 
by the compressive perturbations. CR pressure is usually 
neglected in the dispersion relation for the compressive per- 
turbations. However we know that it is there as CRs and 
the perturbation are connected by the magnetic field. When 
a fraction of the CR pressure is lost due to the process de- 
scribed above, some of the perturbation energy is lost too. 
This provides a new mechanism for the damping of com- 
pressive perturbations in a medium with CRs. 

As we noted previously, each particular scale of com- 
pressive motions A transfers energy to a smaller scale of 
motions at Vp. Also, we know that the flow of energy in 
Kolmogorov-type turbulence is a constant. The steepening 
rate for the reasonable range of parameters is, however, 
growing rather fast, as we show in the next section. This 
could lead to a model where compressive motions are fully 
damped at a particular scale Afb where Kolmogorov energy 
transfer rate is equal to the steepening rate at the scale 
'■p(Afb). More careful consideration, however, shows a differ- 
ent picture. At the end of §4 we showed that scales larger 
than A could also provide some compression for the CRs. 
If this compression is higher than the one produced at A, it 
will also provide higher energy rate. Therefore, the compres- 
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Figure 2. Energy density of compressive modes and Alfvenic 
slab- type waves, induced by CRs. The energy is transferred from 
the mean free path scale to the CR Larmor radius scale. If the 
mean free path falls below compressive motions cutoff or feedback 
suppression scale, the spectrum of slab waves becomes steeper. 



sive motions at scales, smaller than our "feedback damping" 
scale are not damped fully, but rather to the state where they 
provide marginal compression. It is easy to see, that this cor- 
respond to the Sv X law, or spectrum of turbulence E~'^. 
This spectrum of compressive motions will protrude from 
the scale Afb down to scale /cut- The amplitudes of the high- 
frequency slab Alfvenic motions, corresponding to these and 
lower scales will be determined by formula (10). The follow- 
ing handy expression and equations from §4 could be used 
to derive Afbi 

(1=) 



Afb 

T 



Fig. 2 shows energy densities of the compressive mo- 
tions and Alfven slab- type waves in two cases: first, when 
the CR feedback is unimportant, but compressive motions 
are damped at Zcut; second, when the CR feedback become 
important at some scale Afb. Spectral slopes correspond to 
the case with fi = 1/3. 



7 ASTROPHYSICAL CONSEQUENCES AND 
DISCUSSION 

7.1 CRs in ISM and galaxy clusters 

Lets consider CRs in galaxy clusters. The magnetic field 
magnitude and the density of CRs are somewhat uncertain 
there (see Esslin et al. 2005), so we adopt values similar 
to our galaxy, namely B — ncR^E > 1 GeV) = 

4 • 10"^" cm~^ and a = 2.6. This corresponds roughly 
to equipartition between CR and magnetic field energies. 
In clusters these energy densities are around 5 per cent 
of the thermal energy density. We will have then Li ~ 
6 ■ 10~^ pc. The reference Larmor radius of 1 GeV pro- 
ton is ro ~ 2 • 10~^ pc. We take the scale L = 1 kpc, 
which, being Alfvenic at this scale, corresponds roughly to 
driving with the virial velocity at the scale of 30 kpc. For 
these numerical values and fi — 1/3 we will have, from 
Eq. JHJ, 5B/B = 0.04(rp/ro)~'''^, almost independent on 
scale, rp,crit ~ lOVo ~ 2 ■ lO"* pc, A ^ W-^{rp/ro)'-'^ pc, 
and the mean free path corresponding to the turbulent 



damping (rp = rp,crit) is 0.4 pc which is much smaller 
than the outer scale. We estimate coUisionless cutoff as 
'cut ~ 10~^ pc. The feedback mean free path will be around 
0.3 pc, so the spectrum of Alfvenic slab motions will be 
mostly steeper, k~^'^ and the mean free paths will be mod- 
ified according to (10) and (5). The efficient CR scattering 
entails efficient second-order Fermi acceleration, see eq. (13), 
the process that may be important for clusters of galaxies 
(Cassano, Brunetti, 2005). 

In our galaxy one can assume same values for B, a 
and n and value of L around 50 pc. We assume an acoustic 
turbulence spectrum for fast waves, taking fi — 1/4. We 
generally get a smaller range of Alfvenic slab motions, from 
scales of ro to about 600ro with SB/B = 0.093(rp/ro)"° ". 
The resulting mean free paths A vary from 2.3 • 10~^ pc 
to 8 ■ lO"'^ pc. In the Galactic Corona, fast waves will be 
damped by coUisionless damping (see, e.g., Ginzburg 1961) 
with a cutoff of around 1.6T0~'^ pc, which is within the range 
of A that we deal with. In the warm ionized medium (WIM) 
the coUisional damping cutoff will be around 10~* pc. The 
feedback mean free path will be around 7 • 10~^. Again the 
spectrum of slab motions becomes steeper and the mean free 
paths of CRs are modified accordingly. 

In § 2 we assumed that the compression factor A is 
larger than va/c. This assumption is satisfied in galaxy 
clusters, as, from the previously adopted values and n « 
10"^ cm-^ va/c « 10"^ while compression factors for 
scales between 2 pc and 2 • 10~^ pc are between 0.12 and 
5 ■ 10"^. For the Milky Way ISM this condition is satisfied 
much better, as for n « 1, va/c ^ 4 ■ 10~^ , and compres- 
sion factors are generally larger, due to the fact that the 
minimum X/L is smaller. 

The slab Alfven modes had been a part of the CR 
paradigm from the very start of the research in the field (see 
Jokipii 1966). Together with anisotropic components they 
are part of some of the modern models of CR propagation 
(see Zank & Matthaeus 1992, Bieber et al. 1994, Shalchi et 
al. 2006). In our model the slab plane-parallel modes emerge 
naturally as the result of the interaction of compressible tur- 
bulence with CR. Although this mechanism is different from 
the earlier considered processes, it may justify some of the 
earlier calculations invoking slab modes. Unlike earlier the- 
ories we predict the dependence on the amount of the slab 
mode energy on the relative pressure of the CRs. 

As we see, in both clusters of galaxies and ionized gas 
in Milky Way the instability within CR fiuid limits the CRs 
mean free path. Like in scenario discussed in YL04, where 
the compressible fast modes were identified as the major 
CRs scattering agent, compressive modes are essential for 
scattering. However in this treatment, unlike YL04, we show 
that compressions at scales much larger than the resonance 
scale are important. This difference is crucial for scatter- 
ing of low-energy CRs, as the fast mode have coUisional or 
coUisionless cut-offs which, depending on the media, may 
be larger than the low-energy CR gyroradius. In this case 
YL04 appealed to Transient Time Damping (TTD) pro- 
cesses, which are less efficient for scattering than gyrores- 
onance^ . Our present work shows that the slab Alfven mode 



^ In fact the gyroresonance instability can be the major source 
of isotropization during the TTD acceleration. 
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discussed in the present paper can be responsible for effi- 
cient scattering. Another important difference from YL04 is 
that our new mechanism require relatively large total pres- 
sure of CRs (see (|1J). In the case when the pressure of CRs 
is negligible the fast modes could stay the major scattering 
agent (see Petrosian et al, 2006). 

Our model predicts rather small mean free paths, but 
this does not contradict the estimates on the average lifetime 
of the CR in the Galaxy. These lifetimes are estimated to 
be around Galaxy thickness divided by the Alfven velocity, 
which is a powerful support for the models with streaming 
instability. Our model will predict similar lifetimes, because 
it includes turbulence which advects CRs on outer scale 
comparable with Galaxy thickness with velocity of around 
Alfven velocity. In fact, the turbulence itself could be gen- 
erated on these scales by Parker instability. 



7.2 Partially Ionized Gas 

Previous discussion is also applicable to partially ionized 
gas, if the degree of ionization is larger than ~ 90%. Indeed, 
for such high ionization degrees the Alfvenic turbulence cas- 
cades to scales less than the ion-neutral decoupling scale (see 
Lithwick & Goldreich 2001). 

If, on the other hand, the degree of ionization is lower, 
we assume that Alfvenic turbulence is fully damped by ion- 
neutral collisions at the scale Zdamp, and it would not be able 
to provide turbulent damping for fcj_Zdamp < 1. As we saw 
in § 5 the damping for slab waves with fcy is provided by 
turbulent eddies with k± ~ '^y^^' therefore, our slab-type 
component arising from CR instability will protrude up to 
scales as large as IL}^^ . This scale could be substan- 

tially larger than the r^^crit derived in § 5. 

According to Lazarian, Vishniac & Cho (2004) the 
regime of viscosity damped turbulence emerges for Alfvenic 
turbulence at scales less than idamp. This regime is char- 
acterized by a shallow fc~^ spectrum of magnetic pertur- 
bations and it persists down to the ion-neutral decoupling 
scale where it reverts to intermittent Alfvenic turbulence 
that involves only ions. The detailed treatment of the inter- 
actions of CRs with turbulence in partially ionized medium 
is beyond the scope of this paper, however. 



T.3 Thermal plasma mean free paths in galaxy 
clusters 

In the paper above we considered the CR component of the 
ISM or ICM, which are the high energy particles that in- 
teract with the rest of the medium via the magnetic fields. 
These particles have a power-law distribution that arises 
from the acceleration terms that are proportional to the CR 
momenta. The astrophysical plasma, on the other hand, is 
assumed to have a Maxwellian distribution and provides us 
with both conductivity and mass density which are required 
for a MHD treatment. In a fully ionized plasma particle- 
to-particle collisions are Coulomb scattering and the rate of 
the collisions becomes smaller with temperature. With high 
temperature and small density these mean free paths can be 
huge. For example, in galaxy clusters it could be as large as 
4kpc. This lead to apparent contradiction, as particles with 



such a huge mean free path will be subjected to acceleration 
and will not be Maxwellian. 

Schekochihin and Cowley (2005) proposed that thermal 
particles will be scattered by instabilities. They considered 
hydrodynamic as well as kinetic instabilities and considered 
the evolution of the cluster from initial state with no mag- 
netic field. Their argument is that the Reynolds number, 
being initially very low, will increase with increasing mag- 
netic field and the dynamo will self-accelerate. They predict 
folded magnetic fields due to high-Prandtl number dynamo 
and their mean free paths are between viscous scale and the 
reversal scale. 

In this subsection we estimate mean free paths of ther- 
mal particles in a way similar to the rest of our paper, keep- 
ing in mind that there are quite a few other plasma effects 
and some MHD dynamo effects that might be important, so 
that these estimates are still rather speculative. This may 
be excused by the fact that thermal mean free path, viscos- 
ity and thermal conductivity are very important for cluster 
dynamics. 

Anisotropic distributions of thermal particles will ex- 
cite waves with inverse wavevectors of the order of thermal 
Larmor radius, tt ~ 10~^ pc as the instability is exponen- 
tially slow for smaller wavevectors (see Mikhailovskii 1975, 
eq. 10.7). All particles will have approximately the same 
mean free path, and the value of 5B/B that provides scat- 
tering will now refer to the total perturbed magnetic field, 
in contrast with its definition in § 4. Apparently the energy- 
transfer arguments of §6 will be most important, as the 
steepening is very fast on thermal Larmor scales. By equat- 
ing steepening and turbulent energy transfer rates we have 
rr/L ^ {SB/Bf, which gives 5B/B ^ 10-^ A « 10"^ pc. 



8 SUMMARY 

All in all, in the paper above we have demonstrated that 

1. Turbulent compressions of magnetic field result in 
the kinetic instability of CRs that drives Alfvenic pertur- 
bations of much higher frequency with wave vectors almost 
parallel to the magnetic field. These Alfvenic perturbations 
efficiently scatter and isotropize CRs. 

2. The above effect is present over the limited energy 
range of the CRs. The high energy cut-off is determined 
by the ambient Alfvenic turbulence. The non-linear back- 
reaction via limiting the CR mean free path and the steep- 
ening of the generated waves control the intensity of the 
new slab-type Alfvenic component. This intensity depends 
on both the amplitude of the compressible perturbations and 
CR pressure. 

3. The presence of linear damping of compressible mo- 
tions or the strong feedback damping effect modifies the in- 
stability and results in a slightly steeper spectrum of gener- 
ated Alfvenic perturbations. 
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APPENDIX A: CYCLOTRON INSTABILITIES OF COSMIC RAYS IN PLASMA 

We follow the standard procedure of deriving the dispersion relations of electromagnetic waves in plasma. The field of the 
wave creates a perturbation /i in the particle distribution function fo. We define the current density of the perturbation as 
ji = CTikEk where aik is the conductivity tensor and 

4m 

Uj = Oij H aik 

is a dielectric tensor. The perturbation eigenmodes are determined by the so-called dispersion equation 



= 



Symmetries of the dielectric tensor are determined by symmetries of the initial particle distribution function. In our 
treatment we consider a two-component medium in which most of the contribution into the dielectric tensor comes from the 
thermal isotropic plasma, while the small contribution from CRs is responsible for the instability. Let us consider transverse, 
circularly polarized waves with wavevector parallel to the magnetic field. The dispersion relation will reduce to en ± iei2 = 
(?}? luP' , and it has a large component e^^ = e^°^ ± i^^^ t^i*^ comes from the contribution of thermal plasma and the small 
component e^^^ that comes from CRs. For hydrodynamic waves with wavelengths much larger than the Larmor radius in 
thermal plasma we have e^?' = (? Iv\, where va is the Alfvenic velocity. Such waves propagate along the magnetic field with 
velocity va- In thermal isotropic plasma both circular polarizations of the wave have the same speed, in other words, linearly 
polarized waves are also eigenmodes in sueh plasma. With the introduction of a CR contribution to the dielectric tensor this 
degeneracy could be broken. For example, for the isotropic distribution of CRs, "shifted" from the origin by some streaming 
velocity, the eigenmodes are linearly polarized waves which are unstable if the streaming velocity is larger than the Alfvenic 
velocity (see, e.g., Kulsrud & Pearce, 1969). If the distribution is not shifted, but oblate or prolate the eigenmodes are the 
circularly polarized waves which are either stable or unstable (Mikhailovskii, 1975, Kulsrud 2004). 

In this paper we study the oblate or prolate distributions of CRs which come from the conservation of adiabatie invariant 
for coUisionless particles. The CR component e^^^ is responsible for a small imaginary part in the solution of the dispersion 
relation for uj. Depending on its sign the instability constitutes the growth or damping. It could be shown that in the limit of 
va/c<^ 1 this growth rate is equal to ^ 

2 2 fvlfdF «|, dF\ 3 

where ujc ~ eB/mcy is a particle gyration frequency and F{p\^,p±) is a distribution function of CRs. The ± sign 
correspond to the two MHD modes. We shall limit ourselves to to the plus sign as explained in § 3. 
We introduce a small anisotropy factor A and the unperturbed distribution function Fo as 

. P±-P\\ „ , , / 2 , 2 \ -a/2-1 



where a is introduced in § 2. We assume a to be between 2 and 3 as for the CR distribution in our galaxy. The oblate 

distribution then will be described as 

F(j>)^{pl+pl{l-Af 

and we can, in the linear order to A, calculate that 

Now the expression for the instability rate will be 

TT^e^ncfl Va {-a - 2) A f ujc v± f k^p^ \ 3 

7CR = ^ / FqS -^'-^ +i^c ] dp, 

m c UJC J ncR c y rn J 

where we replaced F with Fo and introduced the cyclotron frequency uic, and the total density of CRs ncR- The integral 
in this expression is dimensionless. The total density of CRs is mostly determined by the low-energy cutoff of the distribution 
F() and is rather irrelevant for the instability where only resonant particles contribute. It is more useful to introduce the 
number of fast particles ncR(p > mLOB/k\\) which is determined by the integration of Fo over a region with momentum larger 
than the resonant momentum. After taking the integrals we denote the gamma function as F and arrive at 

/, ^ ncR(p>ma;B/fe||) tt^/^ r(a/2) 

7cr(A;||) = Wpi —AQ, where Q = ——{a + 2){a - 



32 ^ 'r{a/2 + 3/2))' 



^ Speaking more quantitatively, the degree of anisotropy should be larger than va/c for instabihty to take place (Kulsrud 2004) . In the 
astrophysical section we show that in the typical setting of the ISM or ICM our approximation is accurate enough (see § 6). 
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